Propagation properties and limitations on the attainable entanglement in a driven 

harmonic chain 
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The limitations on the production and profitability of entanglement in a harmonic chain under 
strong driving are considered. We report on the limits of attainable entanglement for a given set 
of squeezings of the eigenmodes, showing that the higher the entanglement the more oscillatory 
and thus less easy to profit from. We also comment on propagation properties of entanglement, 
discussing the role of fast rotating terms and illustrating several issues with the example of a sudden 
switch of the coupling. 
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I. INTRODUCTION 

Entanglement is recognized not only as a striking fea- 
ture of quantum mechanics, but also as an important 
resource for quantum information tasks. Its creation and 
manipulation will be a must if we arc to profit from the 
quantum speed up of some information protocols that are 
not available in the classical world. Many body systems, 
studied under the light of quantum information, have 
provided [l| new tools for the condensed matter commu- 
nity, while they can be regarded as interesting devices for 
the production and manipulation of entanglement. 

A particular example of such systems is the harmonic 
chain, whose entanglement properties have been exten- 
sively studied [ll-Q. In addition to the static proper- 
ties of entanglement in the chain, its production and 
manipulation starting from a ground or thermal state 
have been analyzed in where it was shown that 

parametric changes in the coupling lead to production 
of long distance entanglement. This production can be 
optimized [|| and is robust against realistic dissipation 
from hot environments 0- Furher, the harmonic chain 
can be mapped into experiments with nanomechanical 
resonators [8| and ion crystals in multitrap arrangements 
[i| , and has the exceptional features of being exactly solv- 
able and having an exact measure of entanglement [Toj . 

In this paper we present a detailed analysis of the en- 
tanglement properties of a harmonic chain with nearest 
neighbor time dependent coupling and give expressions 
for the maximum attainable entanglement between dis- 
tant oscillators when the chain is initially in the ground 
state. These expressions are obtained for weak and 
strong coupling regimes, and for any amount of squeez- 
ing in the eigenmodes. It is shown that a high amount 
of squeezing can lead to higher entanglement but at the 
expense of it being strongly oscillatory. This is explained 
in terms of the optimality of angular relations between 
different eigenmodes, and how they affect the usefulness 
of the squeezings to produce entanglement. We show how 
these optimal relations shrink to smaller sets for higher 
squeezing or coupling, and quantify this transition. The 



attainable entanglement is explicitly calculated for the 
simple case of a sudden switch of the coupling. 

The propagation properties in the chain have been 
studied for example in [5[. We show here that reason- 
ing in terms of the group velocity in the continuous limit 
yields correct results for the speed of entanglement, as 
well as insight on the dispersive properties of the chain 
as a medium for excitations. We relate the importance of 
fast rotating terms to the strength of the coupling, as well 
as to the ability to produce entanglement with an opti- 
mal modulation, and compare the harmonic chain to the 
anisotropic XY spin chain where the dispersion relation 
has more parameter freedom. 



II. HARMONIC CHAIN 

The system under consideration is a closed chain of 
./V harmonic oscillators with identical frequency loq and 
harmonic coupling between nearest-neighbor c(t): 
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(we will use m = H = 1 throughtout the paper). To 
diagonalize the Hamiltonian ([T]), wc introduce the normal 
mode coordinates Qi and Pi via, 



1 N 



e 2niln/N Qi, 



i=i 

The Hamiltonian of the chain 
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is then that of a set of independent oscillators with time- 
modulated frequencies, uif(t) = Wq+4c(£) sin 2 (7r//iV) (we 
have used the property, Q-i = Q[, P_; = pj). The linear 
Heisenberg equations of motion for Pi(t) and Qi(t) can 
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readily be written down as: 

Qi(t) = Q t (0)Y l (t)+P}{0)X l (jt) (4) 

Pi(t) = j f Qi(t) = Q t (0)Yi{t) + P/(0)Xi(t) (5) 

where the functions Xi(t) and Yi(t) are solutions to the 
classical equation of motion x = —uji(t) 2 x with initial 
conditions Yj(0) = 1, Y(0) = 0, X t (0) = 0, Xi(0) = 
1. We will restrict ourselves here to Gaussian states, 
whose complete information is contained in the symmet- 
ric 2N x 2N covariance matrix T. Its elements are de- 
fined as: T qnqm = 2Re(q„q m ), F qnPm = 2Re(q„p m ) and 
r pnPm = 2Rc(p n p m ) . The first moments are of no rel- 
evance to the entanglement properties and we can drop 
their description from now on. If we assume that the 
chain is prepared in the ground state of the nonintcract- 
ing Hamiltonian, the time-evolved elements of the covari- 
ance matrix can be written in a simple way: 
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where we have used wor g „ gm (0) = T PnPm (0)/ujo — <5 Iijm , 
r 9 „ Pm (0) = and (Q r Qt) = {QsQt)S r , s - The latter 
equality holds only in the case where the eigenoscillators 
were uncorrelated at the beginning, and because they 
never interact, will always be. Similar expressions are 
obtained for (q n Pm) and (p n Pm)- The problem has been 
reduced to finding the second moments of a set of inde- 
pendent oscillators with modulated frequency. 

It is known that a time-dependent oscillator is 
squeezed when its frequency is changed nonadiabatically, 
while an adiabatic transformation leads to no squeez- 
ing The frequency modulation effects a Bogoliubov 
transformation in the ladder operators of each oscillator: 
a — > fid + v<v which keeps the canonical commutation 
: 1. They can be parameterized as 
e l2e sinhr, thus yielding 



relations \p,\ 2 — \v 
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for the quadratures of the eigenmodes. The time depen- 
dence of the eigenfrequencies wj, the squeezing parame- 
ters ri and the squeezing angles 9i are all controlled by 
the linear coupling coefficient c(t). 



III. ATTAINABLE ENTANGLEMENT 



Here, v_ is the smallest symplectic eigenvalues of the 
reduced and partially transposed covariance matrix of 
the two oscillators. 

We will consider in this study only opposite oscillators, 
m = N/2 + n, but any other pair could as well be eas- 
ily described with our approach. Opposite pairs have the 
largest possible distance in the chain, and it can be shown 
that they exhibit the largest value of entanglement, due 
to interference between signals coming through the two 
branches of the chain. In addition, their exponent in the 
quadratures becomes e 2m K n - m )/ N = (—1)', thus distin- 
guishing even and odd normal modes. Due to the trans- 
lational invariance, the logarithmic dependence of pairs 
of oscillators doesn't depend on their position n but only 
on their relative distance n — m. 

The symplectic eigenvalues of the 4x4 reduced covari- 
ance matrix T n ^ m can be written explicitly using local 
symplectic invariants fl2j . which remain unchanged by 
operations on only one of the oscillators. A reorganiza- 
tion of the invariants for opposite oscillators reveals that 
the symplectic eigenvalues can be expressed as sums of a 
unique quantity Ai tTn : 
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From eq. © it is clear that the entanglement between op- 
posite oscillators will be highest when Ay « x 2 . This 
implies minimization of A Q dd,odd and A evenje ven and max- 
imization of A dd,cven or A cvcn! odd for all eigenmodes. 
Clearly, this is almost intractable and a numerical anal- 
ysis is here required. 

Of course the dynamical evolution given by the time 
function c(t) dictates the sets {ri(t)} and {8i(t)}. How- 
ever, when modulation is over, and a finite constant cou- 
pling c is kept, the squeezings do not change, but every 
angle oscillates with its own frequency uji . We study next 
which combinations of angles are optimal, in terms of 
production of entanglement betwen opposite oscillators, 
for a given set of generated squeezings. 



The logarithmic negativity [Toj gives an exact quantifi- 
cation of the bipartite entanglement between oscillators 
n and m, 



E N = max(0, - log 2 (|i/_|). 



(8) 



A. Optimal angles at moderate squeezings and 
small coupling 

The eigenfrequencies of the normal modes are symmet- 
ric with respect to N/2. That means that Lu^-i = Wj, 
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lun/2 is the greatest, and lo^ = which we set to 1. 
With that in mind it is clear that we only need to ana- 
lyze the behavior of half the number of eigenmodes. Since 
opposite oscillators have a separated odd-even structure 
in the covariance matrix the first nontrivial case for a 
closed chain is having four oscillators. In that case the 
number of eigenmodes, and of angles, to be analyzed is 
2. 

The regime c << u!q and moderate squeezing was stud- 
ied in [(|; in this regime uj m /uj[ ~ 1. There it was shown 
that there exists optimal relations for the angles: 



restricted to a smaller measure. This happens for the 



= 2fc?r/2 

= (2fc + l)7r/2 



I + m even 
I + m odd 



(13) 
(14) 



with k G Z. 

We have repeated such analysis for different number of 
oscillators, and obtained the same result. Hence maxi- 
mum entanglement is achieved when both even and odd 
oscillators synchronously achieve those relations for the 
angles. The limit on attainable logarithmic negativity is 
thus: 



rpmax 
N 



= ~2 l0g2 



I, odd 




(15) 

Entanglement for opposite oscillators is only nonzero for 
nonvanishing squcezings. 

The optimal angular relation is modified (and thus the 
validity of cq. (|15|) ) when the coupling or the squcezings 
increase. A self-consistency argument (see Appendix [A)) 
shows that the optimal relations for the angles is valid 
up to: 
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In next section we will show that for higher squcezings 
and coupling the optimal angles shrink to a smaller set of 
values. This immediately means that for a set of strong 
squeezings {r s } it will much harder to obtain dynamically 
(by independent rotation of each eigenmode at frequency 
lo s ) the optimal combination of angles, and thus entan- 
glement will be higher but much more oscillatory. This 
seems to be a fundamental limitation of the harmonic 
chain, hard to avoid unless the coupling is switched off 
when the angles reach their optimal values. 

We have checked this validity condition up to a chain 
of eight oscillators for different combinations {77 }, see 
for example figure [TJ Indeed, only when we approach 
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e 2r ' ~>J2i e 2r ' they cease to be valid. 



B. Optimal angles at high squeezings and strong 
coupling 

We have checked numerically that by increasing the 
quantity J2i e2r ' the optimal angular relations get 



opposite regime, i.e. -^p J2i e2n ^ J2i e 2n ■ The angu- 
lar relations (Till)) . (Tl"4"|) still hold, but they are restricted 
to areas near the values 6 n = mr/2, with n oc id even and 
n-even odd. The evolution from one regime to the other 
can be seen in a naive representation in figure Q] (naive in 
the sense that we are equating all angles with the same 
parity and all squcezings). 

In the case of small coupling c/w 2 <C 1, but high 
squeezing, we can obtain the maximum attainable log- 
arithmic negativity by assuming the new optimal angles. 
We find the same maximum entanglement E^ ax as in 
Eq. (|15[) . This means that this expression is exact for any 
amount of squeezing. This extends the validity of eq. 
(|15l) to any degree of squeezing as long as the coupling 
coefficient c is small. 

Though, when the coupling is higher and c/ujq cannot 
be neglected, we can still obtain a perturbative expres- 
sion: 
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where lO, mE means a sum running over I odd and m 
even and viceversa. This equation has been obtained 
from the expression of the local invariants in terms of 
Ai, m when expanded in orders of x = E/ e_2 ' V^ as a 
small quantity. The correction factor a is smaller than 1 
so it diminishes the argument in the logarithm, hence in- 
creasing the entanglement. Thus, eq. (|17|) with a ~ 1 is a 
lower bound of the maximum entanglement which can be 
achieved. The higher order corrections increase slightly 
the predicted maximum achievable entanglement. Also 
notice that approximately we have 7 = 0(x 4 ) (it can 
be seen by taking similar squeezings n ~ R, so that 
\ ~ e~ 2R ) already, a small correction indeed. That this 
expression for the maximum entanglement is so similar to 
the former one shouldn't be a surprise since it comprises 
the best possible combination of exponents. Any other 
combination of angles would yield positive exponents to- 
gether with the negative ones. 

From the known expression for the irreversible work 
dissipated into each eigenmode [lH, Wdi ss .,; = sinh 2 77, 
we can conclude Q that obtaining a high level of nonsep- 
arability requires an exponential investment in energetic 
resources. 
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FIG. 1. Logarithmic negativity versus Boddj^ (abscissa) and 
9even/Tr (ordinate) for equal squeezing on all eigenmodes, at 
c = 0.3u>o for eight oscillators. The white peaks mean highest 
entanglement. We observe the transition between the two 
different regimes of optimal angles as squeezing is increased, 

with the quantities ^ E; e 2r ' = (0.09,0.2,0.55,1.5) and 



jv 



-2n 



(0.82,0.37,0.14,0.05). 



C. Transition between regimes 



Notice that since every eigenmode's angle evolves with 
its own eigenfrequency, the evolution of entanglement will 
run over the angular space represented in figure [TJ When 
the area of highest entanglement gets smaller, the time 
spent in this region will be shorter, thus producing a 
more oscillatory character in the time evolution of en- 
tanglement. Concluding, entanglement can be made to 
increase, but it will necessarily oscillate more strongly. 



D. Sudden switch 

The sudden switch is a great example in order to fully 
analyze the attainable entanglement, since we know ex- 
actly the amount of squeezing it produces. An instan- 
taneous change of the coupling from to c, changes the 

eigenfrequencies from ujq to = \J ^ + 4c sin 2 (nl/N). 
Such jump produces a known amount of squeezing, which 
according to our definitions of the variances is 
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If we use this equality to express the logarithmic negativ- 
ity between opposite oscillators, we realize that we have 
reduced the set of parameter to {9i} and c. Eq. (fT7|) 
yields the maximum logarithmic negativity which can be 
achieved for opposite oscillators: 
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that for e.g. eight oscillators yields 



In figH] we show the transition between regimes. We 



see that in the upper left figure, which has e 2n 

0.09 <C J2i e~ 2n = 0.82, the optimal angular relations of 
the first regime hold. Note however that as we increase 
the squeezing, the entanglement begins to be slightly 
higher in regions near the new optimal angles, and finally 
gets highest in regions quite close to the optimal values. 
The second regime of optimal angles doesn't come until 
e2n i s higher than J2i e ~ 2r ' (lower-right figure). 
A more thorough analysis for different combinations of 
squeezings {n}, and independent values for the angles 
confirms the new restricted angular relations which are 
centered around the optimal values. This can be seen in 
the example of figure [2] There we have run over the whole 
space of angles for an arbitrary set of squeezings (fixed 
around a typical value of a given size) and only those 
points with highest entanglement have been drawn (up 
to a deviation of 2%). In the first row we see the regime 
of small squeezing, and in the second raw we show the 
opposite regime. It is clearly seen that the new angu- 
lar relations are the same as before, but restricted to a 
smaller area, as explained before. The same occurs if we 
increase c/w 2 instead of/and the squeezings. 
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This expression yields E^ ax - ~ y-^- for small c/w 2 and 
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\ log 2 c/o; 2 for big c/uiq. Clearly, by increasing 
the value of the coupling constant, there is no limit to 
the entanglement we can produce between opposite os- 
cillators, though it must be noted that the energy cost 
grows linear with c/loq. By using an optimized ramp (see 
e -g- 0) we can though produce very much entanglement 
without the need to reach the strong coupling regime. 



IV. ENTANGLEMENT PROPAGATION 

In Ref. @ the optimal coupling modulation was ob- 
tained for a chain of 8 oscillators, showing maximal en- 
tanglement production between all pairs of opposite os- 
cillators. In this section we would like to remark the 
connection between this optimization and the propaga- 
tion properties of the harmonic chain. In particular we 
will sec that the optimal modulation produces a propaga- 
tion of entanglement along the chain without entangling 
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FIG. 2. Set of optimal angular relationships (in units of n) for c = 0.05o;o and N = 8 oscillators. We have run over all 
angles for a set of arbitrary, but similar, squeezings and drawn the angles with maximum entanglement (up to a deviation of 
2%). In the upper row the set of squeezings is (n = T7,V2 = re,rs — rs,r4) — (0.98,1.07,0.89,0.72) and in the lower row 
(1.92, 2.26, 1.90, 2.37). The same behavior as in in figure [JJ can be seen, tough in this figure we can see explicitly the optimal 
relations between angles of different and equal parity. The lower row shows that for high squeezing the region for optimal 



angles shrinks around the optimal values. The self-consistency condition yields e ' = 0.6(6) and ^ ( 1 

for upper(lower) rows, showing that the condition correctly predicts the change between regimes. 
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oscillators, except for opposite ones. This is in stark con- 
trast with the case of a sudden switch of the coupling Q , 
where the entanglement 'wave' sequentially entangles all 
pairs during its propagation. 



A. Non modulated chain 

We consider here a similar argument as given by Cubitt 
and Cirac fl4j . which states that entanglement can be 
seen as correlations traveling along the chain, the prop- 
erties of which can be controlled through its dispersion 
relation. Such dispersion relation is given by the spec- 
trum of the chain, as we show next. 

Consider the chain in its uncoupled ground state, and 
a sudden switch of the coupling at t = 0. Taking the 
thermodynamic limit N — > 00 and 2-Kk/N —¥ 0, we have 
the continuous spectrum of eigenfrcquencies w(0) = [u)q + 
2c(l — cos^)] 1 / 2 . If we take for example the correlation 
(<7n9m) as representative, we obtain from eqs. ([H]),© and 
<20b: 



of the medium. The group velocity: 



(q n q m ) = 



with x 
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s=±l 



2jv 



d(j)A{(j)) cos [cj>x+2suj {4>)t} (23) 



- m ^ 0, A(4>) = (wg - u} 2 ((t)))/uj 2 ((t>) and 
we have removed a constant term. This is clearly the 
sum of two countcrpropagating wave packets, with dis- 
persion relation given by the spectrum u(<fr). Thus, the 
excitations created by the sudden switch will travel along 
the chain following the dispersion characteristics of the 
medium, given by u>(cj>). From this relation, the group ve- 
locity yields information on the dispersion characteristics 



duj(4>) 



csrncp 



(24) 



is only rather flat when c/wg — > 0. Thus, each wave with 
different momentum will travel with a different velocity, 
causing dispersion of the wave packet. Though this pic- 
ture might seem too naive, it is rather powerful. It can 
be shown that in a more general situation where the cou- 
pling has been strongly modulated during t < 0, but is 
kept constant afterwards (t > 0), the propagation of wave 
packets is also ruled by the dispersion relation uj (</>). 

In Ref. [3l they were able to engineer the parame- 
ters so as to have nondispcrsive wave packets, even in 
a medium with high loss/gain (nonncgligiblc anisotropy 
7). It was so because they had an extra degree of free- 
dom, they could manipulate the coupling strength and 
anisotropy. In our case, the interaction Hamiltonian writ- 
ten in terms of (bosonic) excitations is 



XiXi + \ oc a,ia\ +1 + a^ai+i + h.c. 



(25) 



whereas in their case, the equivalent interaction Hamil- 
tonian is 

(1 + l)°f°i+i + (1 - 7)0? °f+i « a i4+i + JdiOi+i + h.c. 

(26) 

This extra parameter 7 allows for the dispersion relation 
to be rather flat in some cases, leading to nondispersive 
propagation. In our case though, nondispersion occurs 
only when c/uj 2 —> 0, which would mean no propagation 
at all. For finite but small coupling, fast rotating terms 
in interaction picture can be neglected: 
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FIG. 3. Time (abscissa) evolution of entanglement between 
oscillators separated by a distance n — m (ordinate) in a chain 
of N = 40 oscillators after a sudden switch of the coupling. 
Due to translation invariance, entanglement is only a function 
of the distance between oscillators. The contour plot ranges 
from Em = (dark) and En = 0.1 (light) . We can see an 
increase of propagation velocities with higher c, but also an 
increase in the signal's loss. 

so there is no loss/gain. However propagation will be dis- 
persive unless an interaction of the form (1 + + 
(1 — "f)piPi+i is used. In figure [3] we see the propaga- 
tion of entanglement after a sudden switch of the cou- 
pling. For small coupling the packet has no loss, but 
is highly dispersive, while higher couplings increase loss 
during propagation. 

B. Modulated chain 

The situation changes drastically when propagation 
occurs during modulation of the coupling. We would 
like to emphasize that the optimal modulation in (see 
fig. [4]) is very similar to the last figure in Fig. [3] (highest 
coupling), in the sense that the entanglement wave seems 
to be completely destroyed along propagation. However, 
the wave revives completely when it reaches the end of 
the chain, so only entanglement between opposite oscilla- 
tors is achieved. The optimization procedure surely has 
to do with improving the time evolution (in interaction 
picture) of the loss/gain terms, which are also responsible 
for the creation of squeezing, such that the entanglement 
is created and delivered only to oscillators sitting in op- 
posite positions. This exclusivity comes from the fact 
that we have maximized the functional in eq. (|15l) . de- 
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FIG. 4. Propagation of entanglement for an optimal modula- 
tion of the coupling (see Appendix [B]| in a chain of N = 8 os- 
cillators as a function of time. The coupling has been switched 
from c = to c = 0.5uq. The quantity n — m is the distance 
between oscillators. Only opposite oscillators (and slightly 
also the immediate neighbors) get entangled, i.e. the wave 
seems to get dispersed but manages to arrive to the end in its 
full form. 



signed to be the maximum achievable entanglement for 
opposite oscillators. Had we wanted to obtain maximal 
entanglement between any other pair, we would have got- 
ten a different expression, with a different combination of 
exponents. Maximizing that other expression we would 
have come to a propagation in which entanglement is only 
delivered to the wanted pairs, if such an optimal modu- 
lation exists at all (remember that the case of opposite 
oscillators is highly symmetric, and it might be the case 
that we cannot deliver entanglement to nonopposite os- 
cillators without entangling some other pairs) . Basically, 
changing the functional to optimize equates to changing 
the weight given to the different 77 , thus favoring delivery 
(once modulation has stopped and the coupling is kept 
constant) of the 'sum of squeezings' to pairs of oscillators 
with a given selected distance n — m. 

The role of the loss/gain terms, which destroy /create 
excitations, in the production of distant entanglement is 
clearly highlighted e.g. for a driven anisotropic XY spin 
chain [TBj . where a resonant modulation of the coupling 
is equivalent to the limit 7 — > 00 in (|26[) . That is, it 
maximizes the presence of creation terms in interaction 
picture. 



C. Transmission speed 

Finally let us comment on the validity of the group ve- 
locity eq. (f2~4~)) . If we consider the simple picture in which 
the highest group velocity is the one responsible for dis- 
tant oscillators to start becoming entangled, we can re- 
strict the analysis to 

tW = «(^ = 7r/2) = -=^=. (28) 
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FIG. 5. Entanglement between opposite oscillators when the 
coupling is suddenly switched from to c = 0.05ojo (black), 
c = O.Iwq (gray) and c = 0.2o>o (light gray). A speed up of the 
synchronization process by a factor ~ 3.5 is observed. 

So propagation can be sped up just by increasing the 
coupling. In figure [5] we show the case of a sudden switch 
in the coupling, where c/wq has been increased from 0.05 
to 0.2; the arrival time of the first peak of entanglement is 
a factor ~ 3.5 smaller for c = Q.2ujq than for c = 0.05wq. 
Using our expression for the group velocity we obtain a 
factor 3.54 . On the other hand, according to the detailed 
study [j| the time at which two oscillators at a distance n 
are not anymore separable is r ~ ny/ojQ + 2c/2c, where 
n is the distance between oscillators. Their expression 
coincides with our very simple expression coming from 
the group velocity r = n/w max . 

Hence we see that reasoning in terms of a dispersive 
medium for excitations, correlations and entanglement 
yield quite similar insights. Furthermore, it allows to 
make quite accurate predictions from very simple argu- 
ments. 



V. CONCLUSION 

We have considered the limitations on the production 
of entanglement in a driven harmonic chain for differ- 
ent regimes, deriving expressions for the maximum at- 
tainable entanglement between opposite oscillators given 
a set of produced squeezings in the eigenmodes. The 
validity of those expressions have been investigated for 
strong/weak coupling and moderate/high squeezing. 

We have also shown that the optimal phase relations 
for the eigenmodes, which provide maximum entangle- 
ment for a given set of squeezings, shrink when the 
squeezings increase, leading to a highly oscillatory be- 
havior of entanglement in time. This imposes a practi- 
cal limit on how much entanglement can be used in a 
harmonic chain, unless we are able to freeze its evolu- 
tion (switching off the coupling) at its maximum value. 
The validity of the phase relations for weak coupling and 
moderate squeezing has been estimated (see Appendix 
|A"|) . We have exemplified the transition to the smaller set 
of optimal phase relations in figs. Q]and[5] 



Finally, we have investigated the role of loss / gain terms 
in the interaction Hamiltonian, relating their effect to the 
validity of the RWA and the strength of the coupling. 
These terms are negligible for weak coupling, leading to 
lossless (though dispersive) propagation. It is precisely 
these terms, which also create excitations in the chain, 
that are favored by an optimization of a time dependent 
coupling. This creates squeezing and restricts delivery 
of entanglement to the selected pair of oscillators. The 
concept of a chain as a dispersive medium for correla- 
tions has proven fruitful for understanding the propaga- 
tion characteristics, which we have used to derive in a 
simple fashion the time it takes for distant oscillators to 
become entangled. 
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nich (NIM), and partially by FISICOS(FIS2007-60327), 
CoQuSys(200450E566) and DiSQuC(AAEE0113/09) 
projects. We thankfully acknowledge R. Zambrini for 
manuscript revision. 



Appendix A: Validity regime for optimal angles 

We will follow here a self-consistency approach in order 
to deduce the regime of validity of the optimal angles 
given in section III. We saw that for low squeezing the 
full relations: 

0i-6 m = 2kTr/2 , Z + meven (Al) 

9i - 9 m = (2k + l)7r/2 , Z + modd (A2) 

with k € Z, are optimal, while for high squeezing the 
range of optimal angles is reduced to the smaller set: 

9 n = mr/2 with n dd even and n even odd (A3) 

The first relations will not hold if the squeezing and/or 
coupling is increased. The angular relations which do 
not overlap with regions near the angles in eq. (|A3|) . 
for example the angles 9 S = 7r/4(mod 71") (s even) and 
6 S = 37r/4(mod tt) (s odd), will be the first to become 
nonoptimal (see fig. [T]). Our strategy will be to compare 
the symplectic eigenvalue in © for the latter angles and 
see when they give a higher eigenvalue than the one given 
by the actually optimal angular values in (|A3[) (i.e. a 
lower attainable entanglement). For the first regime we 
choose 9 even = K / A and 9 dd = 37r/4, while for the other 
regime we choose 9 even = n/2 and 9 a dd — (note that 
this retains the generality of the argument). Following 
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the notation of Eq. (j9)) we can write: 



,2(r ( +r m ) 



lO,mE 



(A4) 




and 



,(f .0) 



(A5) 



(l+^Eu^e 2 ^" (A6) 



=(l + f)(E i0 ^)(E m0 e 2 -)x (AT) 



(E iB e 2 -)(E 



(A8) 



where we have approximated ^ ~ (1 + a/,m^) (that 

is, c/tu 2 <C 1) and taken a; )TO of order 1. We use also the 
previous notation where lO, mE mean I over odd values 
and m running over even values. 

The mixture of sums for odd, even values of these terms 
makes it difficult to obtain a clear picture. Therefore 
we employ a further assumption, that the squeezings are 
similar for even-odd cigenmodes in the sense that if we 
define 



B 



we approximate 



O t (E) 



E 

lO{E) 



,±2r, 



(A9) 



(A10) 



This is reasonable in the sense that c(t) is the only con- 
trollable parameter in the system, which squeezes all 
cigenmodes even if it is tailored to only squeeze a given 
parity. We further notice that B~/B + <C 1 whenever 
the squeezings are relative strong (e.g. r > 1/2). 

With these assumptions we can now compare the sym- 
plectic eigenvalues of each regime to conclude that 



\S^ ] \>\^ fi) \ 



, i.e. when 



E^>E< 



B+ > B~ 



-2rj 



(All) 



(A12) 



We have checked, up to a chain of eight oscillators, 
that the optimal angles are valid, according to the above 
relation, for different combinations {r s }. They cease to 
be valid only when we approach j% J2i e2r ' ^ Ez e 



-2r, 



Appendix B: Optimal control algorithm 

Optimal control theory (OCT) is a branch of applied 
mathematics whose objective is the determination of the 



function which optimizes a given cost functional (l6j ] - A 
very famous particular case leads to the Euler-Lagrange 
equations of motion in classical mechanics, when the ac- 
tion functional is extremized with respect to trajectories. 
We derive the OCT algorithm which has been used here 
and in [(| [l3| , for the sake of reproducibility. 

In our particular case, we want to find the coupling 
modulation c(t) which maximizes equation (|17|) . In this 
problem, this can be achieved with purely classical equa- 
tions. The solution to the time evolution of the chain is 
fully determined by the equations of motion of the normal 
modes, which are quantum harmonic oscillators with fre- 
quency modulation. Examining eqs. ([7]) we observe that 
the average energy of a given normal mode is 



■ cosh2r. 



(Bl) 



therefore the squeezing of a given mode can be ex- 
pressed purely by means of the quantity 2{E s )/u s . or 
in terms of (p 2 ) and (x 2 ). Thus, for a given modulation 
of the coupling coefficient c(t) the squeezing of eigen- 
mode s can be obtained from the classical equation of mo- 
tion of a frequency modulated oscillator with frequency 

^ to 2 + 4c(i) sm 2 (Trs/N) and arbitrary initial conditions, 
provided that its initial energy is equal to the ground 
state energy of the quantum oscillator. 

The aim of the optimal control algorithm is to min- 
imize the cost functional J (c(t)) of the coupling coeffi- 
cient in equation (fT5|) . In our case this functional de- 
pends only on the coordinates of the modulated clas- 
sical oscillator at the total modulation time r, hence 
J(c(t)) = h(x(T)). The vector x = (afi, ...,xn) con- 
tains the positions and momenta of the eigenmodes such 
that Xi = (xi,pi). A simple way to impose the equa- 
tions of motion of the eigenmodes is to include Lagrange 
multipliers in the cost functional, so at the end we have 



h(x(T)) 



dt 



dh. 
dx 



a(x, c) 



(B2) 



where £ = (£i, £tv) are the multipliers (with £j = 
(xf,pf) ), or costates in optimal control jargon, and 
the equations of motion of the eigenmodes are simply 
x = a(x, c) . In order to find the equations governing the 
optimal control algorithm, we need to follow an approach 
very similar to the deduction of the Euler-Lagrange equa- 
tions, which can be found in[l6|. The resulting equations 
are: 



1) 

2) 

3) 



IV) 



Si (a?, c) 
0(1*. c 



dxi 



dh 
dx 



4) £* -a(x*,c*) <C -a{x*,c) , Vc 

where 1) is the evolution equation of the eigenmodes, 2) 
is the evolution equation of the costates, 3) are the initial 
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conditions of the costates (note that they are stated at 
the final time) and 4) is Pontryagin's minimum princi- 
ple which states that the optimal trajectories are those 
which minimize £ ■ a all along the trajectory and the stars 
mean "optimized" (it should be noted here that Pon- 
tryagin's minimum principle is an extension of optimal 
control theory, typically used when the control param- 
eter c(t) is bound in a given range). The fact that the 
costates have "initial conditions" at the final time pre- 
vents us from obtaining the optimal c(t) in one go. There 
are several alternatives for solving the problem, but we 
have chosen to use an iterative scheme with steepest de- 
scent as follows: 

a) Choose a trial function c(t). 

b) Evolve the eigenmodes until r (and record their tra- 
jectory). 

c) Obtain, through equations 3), the initial conditions for 
the costates. 

d) Evolve backwards in time the costate through their 
equations 2) (and record their trajectory). 

d) Change the old c(t) by the amount a<9(£ • a)/dc. 

e) Repeat the process starting with the improved c(t). 



The factor a has to be changed according to the perfor- 
mance of the algorithm. Too big an a will tend to find 
too fast the solution, which we will recognize as an oscil- 
lating solution. If a is too small, the convergence will be 
slow. 

In our system the equations of motion for the eigen- 
modes are ii = pi, pi = —ujfxi, while the equations 
for the costates arc xf = uifpf and jq = —x\. Wc 
also have a = (p, —uj 2 x) and thus the gradient function 

<9(£ • a)/dc = -4 J2k sil>2 (w) Xi Pi- Thc cost mnct i ona l 
can be chosen to be the argument in the logarithm in 
cq. (fT7|) . the inverse of the sum of energies, etc. A bit of 
trial and error is needed until a suitable cost functional 
is found. 

A special comment should be made here. It must be 
noticed that the role of the cost functional h in the im- 
provement in c(t) is quite hidden in the algorithm. Its 
effect is felt only through the initial values for the evolu- 
tion of the costates, which at the end affect the gradient 
function <9(£ • a)/dc. 
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